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ABSTRACT 
Water  quality  data  is  very  important  in  a  water  management 
program.    The  data  is  used  to  detect  changes  in  water  quality,  to 
ascertain  the  amount  of  chemicals  flowing  into  lakes,  to  give 
detailed  information  on  the  cause  and  effect  relationships  of  the 
factors  affecting  water  quality.    Such  survey  programs  are  expen- 
sive to  operate  so  attempts  must  be  made  to  eliminate  the  necessity 
of  taking  some  samples  which  are  producing  similar  water  quality 
information  and  still  have  data  adequately  describing  water  quality 

variations. 

In  this  report,  spatial  and  temporal  factors  affecting  water 
quality  are  discussed.    Analytical  statistical  techniques,  both  para- 
metric and  non-  para  metric,  are  presented  to  determine  whether  there 
are  significant  differences  between  data  of  consecutive  months.    The 
selection  of  parameters  and  groups  is  discussed  and  the  techniques 
are  applied  to  data  from  Duffln  Creek.    The  results  are  presented  and 
various  recommendations  are  made. 
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INTRODUCTION 

With  the  growing  demand  for  good  water,  it  is  necessary  to 
determine  the  quality  of  existing  water  supplies  in  the  different  areas 
covered  in  a  water  management  program  and  to  assess  the  changes, 
if  any,  taking  place  in  the  water  quality.    The  efficiency  and  adequacy 
of  the  monitoring  program  should  be  given  particular  attention  in  order 
to  ensure  a  sufficiently  large  supply  of  data  and  to  indicate  what 
measurements  should  be  made  in  future  monitoring  programs.    To 
carry  out  the  monitoring  program  economically,  with  existing  survey 
techniques,  it  is  necessary  to  determine  how  often  and  how  many 
samples  should  be  taken  for  the  analysis  and  interpretation  of  data. 

The  Water  Quality  Surveys  Branch  is  conducting  a  regular 
water  quality  monitoring  survey  program  of  the  river  basins  in  the 
Province  of  Ontario.    The  purpose  of  the  program  is  to  ascertain  the 
amount  of  various  chemicals  flowing  into  the  Great  Lakes,  and  to 
develop,  based  on  the  survey  data,  a  program  in  the  watersheds  for 
maintaining  clean  water.    Collected  over  several  years,  the  water 
quality  data  can  be  used  to  indicate  long-term  environmental  changes 
in  watersheds  and  to  give  detailed  information  on  the  cause  and  effect 
relationships  of  the  factors  affecting  water  quality.    Effective  water 
management  rests  on  a  sound  water  quality  monitoring  system.    Since 


such  systems  are  expensive  to  operate,  It  is  important  that  the  sam- 
pling program  Is  designed  to  produce  sets  of  data  which  can  be  inter- 
preted using  reliable  statistical  techniques  and  is  the  most  efficient 
one  possible. 

The  question  often  asked  is  how  many  observations  should  be 
made  to  ensure  that  the  risk  of  making  a  wrong  decision  is  acceptably 
small.    Some  risk  is,  of  course,  always  present  because  of  the  random 
variation  in  results,  but  the  risk  becomes  smaller  and  smaller  as  the 
number  of  observations  Increases.    It  is  clear,  however,  that  for 
reasons  of  economy  of  time  and  effort,  the  extent  of  the  monitoring 
program  should  be  kept  to  a  minimum,  consistent  with  the  maximum 
risk  we  are  prepared  to  accept  of  making  a  wrong  decision.    In  practice, 
there  will  be  a  certain  optimum  rate  of  sampling.    However,  it  is  often 
better  to  evaluate  the  cost  of  sampling  and  the  economic  implications 
of  making  a  wrong  decision  before  one  can  consider  decision- making 
on  the  basis  of  a  specified  risk  of  being  wrong.    A  partial  answer  to 
this  question  is  found  in  the  theory  of  sampling  which  will  be  discussed 
here  briefly. 

Before  commencing  on  the  adequacy  of  background  data  or  on 
data  which  should  be  developed  in  the  future,  it  is  necessary  to  evaluate 
the  existing  sampling  data.    This  evaluation  can  help  to  determine 
whether  the  data  of  particular  months  and  years  is  or  is  not  contributing 


I 
I 


significantly  to  the  total  water  quality  Information.    It  can  also  indicate 
whether  or  not  it  is  possible  to  eliminate  the  necessity  of  taking  some 
samples  which  are  producing  similar  water  quality  information  and  still 
have  data  which  adequately  describe  water  quality  variations.    There- 
fore, this  study  was  initiated  to  test  various  known  statistical  tech- 
niques for  their  applicability  in  detecting  similarities  and/or  dif- 
ferences in  the  water  quality  information  collected  during  the  five 
year  period  (1964-68). 
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I    VARIATION  OF  WATER  QUALITY  DATA 

Space  and  time  are  two  elements  which  affect  all  important 
water  quality  constituents.    Therefore,  an  insight  from  this  point 
of  view  may  provide  some  information  on  the  optimal  monitoring 
conditions, 
(a)     Spatial  Variation 

The  quality  of  water  in  a  flowing  stream  varies  from  place 
to  place  within  the  system,  chiefly  as  a  function  of  the  physio- 
graphic and  geological  characteristics  of  the  tributary  drainage, 
the  distribution  and  type  of  man-generated  and  natural  sources  of 
upstream  pollution  and  the  hydrologic  properties  of  the  basin.    In 
general,  it  can  be  said  that  the  quality  of  stream  water  is  best  at 
the  highest  point  of  the  drainage,  where  it  is  closest  in  time  and 
space  to  its  original  form  as  precipitation,  and  has   had  minimum 
exposure  to  man.    As  the  water  moves  downhill  and  the  tributary 
drainage  area  is  expanded,  the  potential  for  pollution  increases 
both  through  natural  causes  and  human  activity,   since  the  number 
and  variety  of  pollutants  would  be  expected  to  Increase.    With  the 
introduction  of  pollutants  of  domestic  and  industrial  origin,  usually 
centered  at  points  within  the  system  rather  than  dispersed  as  in  the 
case  of  natural  pollution  input,  the  stream  experiences  abrupt 
changes  in  its  quality.    The  greatest  rate  of  change  occurs  close  to 


pollution  sources.    These  changes  may  be  attenuated  somewhat 
because  of  mixing  with  greater  quantities  of  dilution  water  from  un- 
polluted sources  as  the  drainage  area  is  expanded.    Generally,  it 
can  be  assumed  that  spatial  variation  in  the  quality  of  streams  is 
one-dimensional  in  most  instances.    Changes  along  the  vertical 
or  transverse  axes  are  relatively  minor  in  comparison  with  changes 
along  the  axis  of  flow, 
(b)     Temporal  Variation 

Temporal  variations  in  stream  quality  are  due  to  a  large 
number  of  time- dependent  factors,  most  of  which  are  related  either 
to  the  annual  hydrologic  cycle  or  to  the  various  monthly,  weekly  or 
diurnal  cycles  possible.     Some  of  these  factors,   such  as  precipitation, 
are  predominantly  random  in  character  and  consequently  exert  a  more 
dramatic  influence  on  quality  than  will  be  observed  with  strictly 
periodic  factors.    Because  of  the  random  nature  of  these  factors, 
related  quality  responses  may  be  most  appropriately  evaluated  in 
statistical  form. 

It  should  be  noted  that  hydrologic  factors  and  quality  res- 
ponses associated  with  them  are  composites  of  random  and  periodic 
components,  the  degree  of  randomness  varying  with  time  and  space 
scales.    As  the  scales  increase,  random  effects  become  less  important 
and  the  responses  to  periodic  factors  tend  to  be  more  significant. 


Thus,  it  can  be  appreciated  that  the  suitability  of  statistical  techniques 
in  a  water  quality  monitoring  program  in  rivers  and  streams  is  not 
simply  a  matter  of  arbitrary  choice.    Careful  consideration  must  be 
given  to  the  size  of  the  sample,  the  time  interval  and  the  nature  of 
hydrologic  phenomena  which  either  directly  or  indirectly  affect 
quality.    One  can  say  that  in  the  headwaters  of  a  stream,  where  the 
degree  of  natural  or  artificial  pollution  is  likely  to  be  minimal,  the 
diurnal  cycle  may  not  be  as  important  as  the  annual  cycle.    Low  in 
the  drainage  area  in  a  biologically  more  active  environment  diurnal 
patterns  such  as  the  daily  dissolved  oxygen  cycle  associated  with 
photosynthesis  may  become  controlling.    Annual  quality  fluctuations, 
closely  tied  in  with  hydrologic  phenomena,  tend  to  be  attenuated 
lower  in  the  drainage  and  increased  in  amplitude  high  in  the  drain- 
age. 


II    ANALYTICAL  TECHNIQUES 

The  statistical  technique  selected  for  this  case  study,  known 
as  analysis  of  variance,  is  probably  the  most  powerful  procedure  In 
the  field  of  experimental  statistics.    The  subject  of  the  analysis  of 
experimental  data  is  so  vast  that  one  can  at  best  touch  upon  some  of 
the  basic  ideas  here. 
(a)      Parametric  Techniques 

In  general,  the  variation  in  sample  values  consists  of  two 
types  of  fluctuations:    rardom  (or  sampling)  fluctuations  and  systematic 
fluctuations.    The  function  of  analysis  of  variance  is  to  find  the 
systematic  variations  between  means  of  different  samples  and  to 
determine  whether  they  are  real  or,  in  statistical  language,  whether 
they  are  significant.    Another  purpose  is  to  determine  how  specific 
sources  of  variation  contribute  to  the  total  variance  of  a  quantity, 
and  to  test  whether  the  effect  of  a  particular  factor  is  real  or  likely 
to  have  arisen  as  a  result  of  random  errors.    In  the  simplest  form  of 
this  analysis,  the  mean  of  two  samples  can  be  tested  for  a  significant 
difference;  that  is,  whether  or  not  both  are  from  the  same  population 
with  a  mean  u.    This  is  done  using  the  t-test  which  can  be  found  in 
any  book  on  statistics.-1 


1 

Freund,  J.E.,   Mathematical  Statistics,  Prentice  Hall,  1962.    p. 266- 269. 


8 

Although  the  basic  computations  for  the  analysis  of  variance 
are  usually  quite  straight- forward  and  described  adequately  in  most 
texts,  the  interpretation  of  the  computations  may  become  difficult  for 
more  complex  situations  since  it  depends  on  the  purpose  of  the 
analysis  and  the  assumptions  made  in  setting  up  the  experimental 
model.    The  underlying  assumption  for  these  techniques  is  that  the 
values  within  a  column  are  normally  distributed.    Generally,  in  the 
analysis  of  variance,  the  null  hypothesis  is  formulated.    This  requires 
a  test  to  ascertain  whether  or  not  the  means  of  the  n  samples  (or 
columns)  within  each  group  are  similar  and  represent  the  same  popu- 
lation.    The  method  used  to  test  this  hypothesis  is  based  on  the 
analysis  of  the  total  variability  of  the  data.    If  the  null  hypothesis 
is  true,  then 

k      n 


k     n  p 


which  is  (nk-1)  times  the  variance  of  all  the  data,  is  due  entirely  to 
chance  variation.    Here  x  is  the  mean  of  all  the  data;  namely, 

k       n 
7      r       ta    I     I    *il 

If  the  null  hypothesis  is  not  true,  then  part  of  this  sum  of  squares 
can  be  attributed  to  differences  among  the  means  of  the  n  populations. 
To  be  more  specific  one  can  write: 

*       n  2  r*£*  2v<  2 

I     I  <*!]-*>  s  I      I  (*IJ  -  V     +n£(x|-x) 

i  =  |     ja|  (=1      j=l  1=1 

(one-way  analysis  of  variance) 

or     S    £>ij-*)2      ■       n£Cx(-K)2     +      kfdj-x)2     +     £    t  (xij  -  x»  -  x,  -  x)2 
i=i   j=i  f  =  i  J=i  i=i   j=i 

(two-way  analysis  of  variance) 


It  is  customary  to  refer  to  the  above  equations  as: 

SST         -       SSE  +  SSB 

SST  =        SSA  +  SSB  +  SSE 

where  SST  is  the  total  sum  of  squares,  SSE  is  error  sum  of  squares, 
SSB  as  sum  of  squares  between  sample  B's,  and  SSA  as  sum  of  squares 
between  A' s.    From  these  equations  one  sees  that  SST,  which  is  a 
measure  of  the  total  variability  of  the  data,  has  been  broken  down  into 
the  sum  of  two  or  three  components.     SSE  measures  chance  variation 
(the  variability  within  the  samples  or  experimental  error)  regardless 
of  whether  the  null  hypothesis  is  true.     SSB  measures  chance  vari- 
ation if  the  null  hypothesis  is  true,  but  it  can  be  attributed  in  part 
to  differences  among  the  population  means  if  the  null  hypothesis  is 
false.    In  the  two-way  analysis  equation,  the  total  variability  of 
the  data  is  broken  down  into  terms  attributed  to  differences  among 
the  A' s  (rows  1  =  1,   2...k)  and  B's  (columns  J  =  l,2...n)  and  the 
chance  error.    In  the  two-way  analysis,  there  are  two  null  hypo- 
theses based  on  the  A' s  and  B's;  that  is,  there  is  no  significant 
difference  among  the  A'  s  and  there  is  no  significant  difference  among 
the  B's.    The  null  hypothesis  concerning  A  is  based  on  the  statistic: 

p        SSA/deorees  of  freedom  between  A' s     -       SSA/(k-l) —    =     (n-l)SSA 
A  "  SSE/degrees  of  freedom  between  B'  s  SSE/(k-  l)(n-  1)  SSE 

and  concerning  Bis  based  on: 
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F       =    SSB/(n-l)  =   (k-l)SSB 

B         SSE/(n-l)(k-l)  SSE 

We  reject  the  null  hypothesis  If  FA  or  FB  are  greater  than  FQ  »(jc_i)/  (n-l)(k-l) 

and  R.    ,  (n-i)     (n-D(k-l)'   (tabulated)  respectively.    In  the  above  model 

we  did  not  allow  for  the  possibility  that  there  might  be  Joint  effects  of 

the  two  variables;    that  is,  so-called  interactions.    The  above  values 

of  F.  and  F     may  also  be  written  as: 

F      =    MSA        >       F  tabulated  (rows)       rejection  of  null  hypothesis 


A 


MSE 


F  MSB       >     p  (columns)    rejection  of  null  hypothesis 

B      '  "MST     '       tabulated 

where  MSA,  MSB,  MSE  are  the  mean  squares  for  A  (row  effects),  for  B 
(column  effects)  and  residual  sources  respectively.    If  the  hypothesis 
that  there  are  no  row  effects  or  no  column  effects  or  both  is  rejected, 
some  other  statistical  method*  can  be  applied  to  get  a  confidence 
interval  on  the  difference  between  two  row  effects  or  two  column 
effects. 

Analytical  techniques  for  comparing  the  means  within  a  group 
and  the  variance  (measure  of  the  variation  from  the  mean)  are  well- 
established  and  documented.    The  techniques  used  in  the  case  study 
are  statistically  strong,  if  an  adequate  number  of  samples  are  avail- 
able within  a  group.    Analysis  of  variance  for  multiple  factors  is 
really  only  an  extension  of  the  two-way  analysis. 


Bowker,  A.  H.  and  Lieberman,  G.  J.,  Engineering  Statistics, 
Prentice  Hall,  New  York.     p. 323  and  295. 
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(b)    Non- parametric  Techniques 

The  parametric  technique  for  testing  whether  several  samples 
have  come  from  Identical  populations  is  the  analysis  of  variance  or 
F  test.    The  assumption  for  the  statistical  model  underlying  the  F 
test  are:    that  the  observations  are  Independently  drawn  from  normally 
distributed  populations;  that  the  populations  all  have  the  same  variance; 
and  that  the  means  in  the  normally  distributed  populations  are  linear 
combinations  of  "effects"  due  to  rows  and  columns;  that  is,  that  the 
effects  are  additive.    However,  if  normality  cannot  be  assumed  and 
if  the  other  assumptions  are  unrealistic  for  the  data  available,  it  is 
necessary  to  investigate  the  non- parametric  techniques  of  analysis. 
This  also  would  increase  the  generality  of  the  findings  by  avoiding 
the  assumptions  made  above. 

Although  there  is  no  underlying  assumption  of  normality  in 
non- parametric  methods,  a  requirement  for  a  sufficiently  large 
number  of  values  exists  in  order  to  obtain  results  comparable  to 
the  parametric  methods.    One  must  also  keep  in  mind  that  If  para- 
metric methods  are  applied  to  non- normal  distributions,  some  invalid 
conclusions  may  result.    Both  the  Mann-Whitney  U  and  Wald- 
Wolfowitz  runs  tests  are  useful  for  a  one-way  analysis  of  two  samples. 
As  mentioned  earlier,  one-way  analysis  is  applied  when  samples  are 
Independent,  that  is  when  there  is  no  row  effect  on  the  values  in  the 
columns. 
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If  one  is  interested  simply  in  testing  whether  two  samples 
come  from  the  same  population  or  whether  they  represent  populations 
which  differ  in  central  tendency  (mean),  the  Mann-Whitney  U-test 
would  be  a  more  powerful  test  than  the  Wald-Wolfowltz  runs  test 
because  it  is  specifically  designed  to  disclose  differences  of  this 
type.    The  runs  test  is  designed  to  disclose  differences  of  any  type 
(central  tendency,  variability,  skewness,  etc.)  and  is  thus  less 
powerful  in  disclosing  the  particular  kind  of  difference  needed  here. 
(This  test  can  be  applied  if  sufficiently  large  samples  are  available). 

The  power  efficiency  of  the  Mann-Whitney  U-test  is  95  per  cent  while 
that  of  the  Wald-Wolfowitz  runs  test  is  not  known.    The  concept  of 
power  efficiency  is  concerned  with  the  amount  of  increase  in  sample 
size  which  is  necessary  to  make  these  tests  as  powerful  as  corres- 
ponding parametric  tests.    Thus,  one  can  avoid  having  to  meet  some 
of  the  assumptions  of  the  more  powerful,  parametric  tests  without 
losing  power  by  simply  choosing  a  different  test  and  drawing  a  larger 
total  number  of  samples  (N).    In  other  words,  by  choosing  another 
statistical  test  with  fewer  assumptions  in  its  model  and  thus  with 
greater  generality  than  theTand  F  tests,    and  by  enlarging  N,  one 
can  avoid  having  to  make  the  assumptions  mentioned  earlier,  and 
still  retain  an  equivalent  power  to  reject  the  null  hypothesis. 
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The  choice  of  non-  para  metric  techniques  is  much  more  limited 
for  the  n  related  or  unrelated  samples  than  for  two  samples.    After 
looking  at  the  data  available  in  each  month,  it  was  decided  that  the 
Friedman  two-way  analysis  for  n  related  samples  (where  n>  2),  and 
the  Kruskal-Wallis  one-way  analysis  were  useful  for  testing  the  null 
hypothesis  that  the  n  samples  have  been  drawn  from  the  same  popu- 
lation . 

Both  tests  are  similar  In  that  ranks  are  used.    However,  the 
Kruskal-Wallis  test  does  not  consider  row  effects  as  the  Friedman 
test  does.    The  tests  assume  that  the  variable  under  study  has  an 
underlying  continuous  distribution.    The  Kruskal-Wallis  method  can 
be  used  where  there  is  a  variable  number  of  samples  in  each  group 
(or  column).    Furthermore,  this  method,  based  upon  comparison  of 
averages,  results  in  a  power  efficiency  of  95.5  per  cent  compared 
to  the  parametric  analysis  of  variance. 
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III    METHOD  AND  DISCUSSION 

The  purpose  of  this  study  was  to  see  if  the  number  of  samples 
taken  each  month  or  year  could  be  reduced,  thereby  achieving  more 
economical  operations.    This  reduction  should  be  made  without  re- 
ducing the  usefulness  of  this  data  in  the  future, 
(a)    Selection  of  Parameters 

The  job  of  selecting  the  reliable  water  quality  parameters 

which  are  important  in  the  monitoring  program  Is  a  difficult  one. 

As  a  general  guide,  each  parameter  must  be  measured  frequently 

and  reliably  and  must  give  a  good  indication  of  pollution  in  the 

streams.    For  these  reasons,  the  following  parameters  were 

selected: 

Dissolved  Oxygen  (DO)  Suspended  Solids 

Conforms  Nitrates  +  Nitrites 
Total  Solids  (NO3  +  N02) 

5-  Day  Biochemical  Oxygen  Total  Kjeldahl  Nitrogen 

Demand  (BOD5)  Total  Phosphates 

The  reliability  and  accuracy  of  these  selected  water  quality 

survey  parameters  is  difficult  to  determine.    A  crude  measure  of  the 

variability  of  the  values  is  the  ratio  of  the  maximum  to  the  minimum 

value.    Alow  ratio,   say,  less  than  15,  would  characterize  a  well- 

r 

regulated  stream  with  comparatively  little  expectation  of  quality 
variation  due  to  hydrologic  influences.    Hence  a  minimal  sampling 
program  may  be  called  for.    In  contrast,  a  ratio  in  excess  of  100 
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would  typify  a  stream  with  large  variations  of  a  random  nature  might 
be  expected.    Such  a  stream  could  require  sampling  as  often  as  once 
a  week  in  order  to  adequately  define  the  random  variables  in  quality 
over  the  annual  cycle.    The  range  and  ratio  are  shown  in  Table  1  and 
Table  2A. 

The  most  difficult  indicator  to  work  with  is  the  coliform  count, 
since  it  normally  varies  by  orders  of  magnitude.    Averages  are  depen- 
dent almost  entirely  on  a  few  high  values  and  the  corresponding  ratios 
are  misleading.    Therefore,  analysis  of  this  data  is  a  little  more 
difficult.    The  other  difficulty  involves  parameters  where  the  minimum 
value  equals  zero.    The  ratio  equals  infinity  and  this  is  unrealistic. 

The  accuracy  of  the  measurements  is  affected  by  low  concen- 
trations (as  in  the  case  of  NO3  +  NOz  and  total  phosphate)  and  the 
limitations  of  the  standard  measurement  techniques  used.    The  tech- 
niques for  all  the  parameters  except  total  phosphate  remained  un- 
changed so  that  the  reliability  and  accuracy  did  not  change.     Methods 
of  measurement  should  be  consistent  so  that  any  possible  type  of 
error  remains  constant.    Whenever  there  is  an  improvement  in  tech- 
nique, further  complications  are  introduced  into  the  analytical  studies. 
Though  general  improvements  were  made  in  the  phosphate  determination, 
it  was  assumed  the  accuracy  of  determining  this  parameter  remained 
the  same. 
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For  this  case  study,  the  water  quality  data  collected  at  Duffin 
Creek  was  chosen  because  of  the  long  records  available.    This  stream 
was  representative  of  both  a  rural  and  an  urban  (industrial)  drainage 
basin.    The  particular  station  chosen  was  near  the  mouth  of  the  stream. 

{b)    Selection  of  Groups 

After  some  thought  it  was  decided  that  the  data  for  each 
selected  parameter  could  be  divided  into  groups  of  months.    Thus, 
months  within  each  group  would  be  compared  with  one  another  to 
see  if  the  data  for  each  month  differed  significantly.    Each  group 
contained  n  months  where  n  could  vary  from  2  ° . .  12,  depending 
upon  the  grouping.    As  the  monitoring  data  was  available  from  1964- 
1969,  each  month  had  at  least  4  or  5  samples.    The  final  grouping 
of  months  is  shown  in  Figure  1.    Each  parameter  usually  had  a 
slightly  different  grouping  pattern  though  general  trends  were 
clearly  visible. 

The  initial  selection  of  the  groups  was  somewhat  intuitive 
and  based  on  trial  and  error.    To  begin  the  selection,  visual  aids 
were  used;    that  is,  plotting  the  monthly  values  and  moving  averages 
for  two  or  three  months.    The  change  in  the  slope  of  the  line  indicated 
some  grouping  possibilities  which  were  used  as  a  starting  point. 

In  checking  and  arranging  the  monthly  data  for  the  two-way 
analysis  of  variance,  it  was  noticed  that  in  the  period,   1964-1969, 
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some  readings  were  missing;  in  other  cases,  there  were  two  or  three 
readings  per  month  (it  was  decided  to  use  only  one  observation  per 
month).    Naturally,  two  questions  arose.    What  values  should  be 
substituted  for  the  missing  values  and  which  value  of  two  or  more 
monthly  values  should  be  taken;  for  example,  the  highest  or  the 
lowest?    For  the  first  question,  it  was  decided  to  substitute  the 
mean  value  for  the  missing  values,   so  that  the  substituted  values 
did  not  greatly  affect  the  comparison  of  group  means,  especially  in 
the  case  of  the  two-way  analysis  of  variance  since  the  two-way 
classification  required  an  equal  number  of  readings  in  each  column. 
As  for  the  second  question,  it  was  decided  to  use  the  highest  of  the 
readings  in  a  particular  month  because  high  values  were  more 
critical  from  a  pollution  point  of  view  in  this  study.    The  one  excep- 
tion was  EX)  where  low  values  are  more  critical  in  pollution  deter- 
minations.   Therefore,  low  values  were  used  in  the  two-way  analysis 
of  variance, 
(c)    Application  of  Parametric  Analysis  of  Variance 

After  the  groups  were  decided  and  the  values  selected,  a  two- 
way  analysis  of  variance  was  applied  to  the  data.    The  two-way 
analysis  of  variance  indicates  the  discrepancies  among  the  means 
of  different  months  and  years  within  each  group;  that  is,  whether 
they  may  reasonably  be  attributed  to  chance  or  indicative  of  differences 


among  the  means  of  the  corresponding  populations.    Alternatively, 
one  can  ask  whether  it  is  possible  to  test  if  two  or  more  sample 
means  (means  of  each  month)  are  likely  to  have  come  from  the  same 
hypothetical  population  with  mean  u  and  variance 

In  this  case  study,  the  months  were  used  as  columns  and  the 
years  as  rows  and  one  reading  per  combination  was  available.    Tables 
1A,  3A,  and  3B  have  been  prepared  in  such  a  way  that  either  of  the 
two  effects,  column  effect  or  row  effect,  can  be  studied  by  two-way 
analysis.    The  important  thing  is  the  degree  of  freedom  shown.    Tables 
1  and  2A  also  show  some  further  computations  including  a  one-way 
analysis  done  on  the  same  data  but  excluding  the  substituted  mean 
values  for  missing  ones  using  an  unequal  sample  method.      This 
information  along  with  the  applied  F  test  was  used  to  show  if  there 
was  any  evidence  to  support  the  hypothesis  that  the  means  of  the 
months  within  each  group  were  significantly  different.    A  significant 
difference  was  indicated  by  SD;  no  significant  difference  by  NSD  in 
all  the  tables. 

Further,  a  one-way  analysis  of  variance  was  also  performed  on 
all  the  data  available,  and  this  required  no  substitution  of  values  since 
the  columns  did  not  have  to  have  the  same  number  of  readings  (i.e. 
using  the  unequal  sample  method).  *    The  results  are  shown  in 

1 

Introduction  to  Probability  and  Statistics,  Alder  and  Roessler, 

Freeman  1969  (Fourth  Ed.)  p.  253. 
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Table  2.    It  should  be  noted  that  a  level  of  significance,  a     ,  of 
0.05  was  used  for  all  the  tests  applied,  both  parametric  and  non- 
parametric  . 

In  view  of  the  lack  of  information  on  the  population  distri- 
bution and  the  unavailability  of  proven  transforms,  it  was  felt  that 
no  transformation  should  be  applied.    Nevertheless,  it  should  be 
borne  in  mind  that  if  the  transformations  were  known,  they  could  be 
applied  to  small  samples  (only  5-10  values)  with  confidence. 

One  further  check  was  made  in  this  study,  though  the  actual 
results  will  not  be  reported  here.    Instead  of  using  the  highest  values, 
the  lowest  of  the  values  for  each  month  were  substituted  Into  the 
final  grouping  for  all  the  parameters  except  DO.    A  full  analysis 
was  done  on  these  data  as  above  for  the  high  values  and  the  results 
were  compared.    It  was  found  that  the  results  were  more  ambiguous 
and  misleading  than  when  high  values  were  used  though  in  most  cases 
there  was  good  agreement.    This  was  particularly  true  where  values 
were  very  low  (especially  nitrate  and  nitrite  (N02  +  N03))in  May, 
June,  July  and  August  where  there  were  a  large  number  of  zero  values 
and  values  of  the  same  magnitude).    This  confirmed  the  original 
decision  to  use  high  values  for  the  two-way  analysis. 

If  a  group  of  three  or  more  months  (columns)  showed  no  signi- 
ficant differences  among  the  columns,  any  combination  of  the  months 
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within  the  group  also  showed  no  significant  difference  for  both  new 
groups, 

(d)    Application  of  Non- parametric  Methods 

Generally,   parametric  tests  are  stronger  than  non- parametric 
tests,  especially  for  given  limited  values.    One  needs  more  values 
to  make  the  non- parametric  tests  as  powerful  as  the  parametric  tests. 
However,  to  test  results  based  on  some  of  the    assumptions  implied 
in  the  use  of  the  parametric  tests,  the  non- para  metric  tests  mentioned 
above  were  used  as  a  check. 

In  applying  the  non- parametric  techniques  to  the  groups, 
several  factors  came  into  consideration.    Two-way  non- parametric 
tests;  for  example,  the  Friedman  two-way  analysis  of  variance  by 
ranks,  could  only  be  applied  In  cases  of  three  or  more  columns  per 
group.    There  was  insufficient  data  for  a  two-way  analysis  of  groups 
with  only  two  columns,   so  that  tests  like  the  Wllcoxon  matched  pairs 
test  could  not  be  used.    For  this  reason,  it  was  discovered  that  the 
one-way  tests  were  more  useful  to  this  study,  especially  since  all 
the  available  data  could  then  be  utilized. 

There  were  three  one-way  non- parametric  tests  used:    the 
Kruskal-Wailis  one-way  analysis  of  variance  by  ranks,  the  Mann- 
Whitney  U-test,  and  the  Wald-Wolfowitz  runs  test.    The  Kruskal- 
Wailis  test  could  be  used  only  for  three  or  more  columns  per  group, 
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but  gave  results  agreeing  well  with  the  corresponding  parametric  tests. 
The  Mann-Whitney  U-test  was  applied  to  groups  with  only  two  columns 
and  the  runs  test  was  applied  wherever  there  was  sufficient  data,  as 
a  check  on  the  Mann- Whitney  test.     Generally,  the  Mann- Whitney 
U-test  was  regarded  as  the  stronger  test. 

In  using  the  non- parametric  tests  it  was  important  to  ensure 
that  the  tests  were  two-tailed.    In  one  case,  for  the  Mann- Whitney 
U-test,  the  table  of  probabilities  (Table  J,  p.271-2732)  given  in  the 
text  was  only  for  a  one-tailed  test.    The  appropriate  values  had  to 
be  doubled  for  a  two- tailed  test.    In  using  any  tables  it  is  important 
to  cneck  whether  they  are  two- tailed  or  one- tailed. 

The  main  difficulty  encountered  in  the  use  of  the  one-way  two 
sample  tests  was  the  occurrence  of  ties;  that  is,  two  or  more  readings 
of  the  same  magnitude.    There  was  no  provision  for  either  deciding  a 
tie  or  considering  it  in  the  calculations.    This  presented  some  dif- 
ficulty in  making  a  decision  on  whether  a  significant  difference 
existed.    It  seemed  best  to  consider  several  possible  arrangements 
but  not  for  all.    After  considering  the  possibilities  it  was  decided 
there  was  actually  no  significant  difference.    In  general,  the  results 
of  the  non- parametric  tests  closely  paralleled  those  of  the  parametric 
tests.    Tables  1A  and  2  shows  the  results  of  all  the  tests. 


2 

Non- parametric  Statistics  for  the  Behavioral  Sciences,  S.  Siegel, 

McGraw-Hill,  p.271-273. 
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IV    SUMMARY  AND  CONCLUSIONS 

In  this  study,  the  most  important  statistical  technique  used 
was  the  analysis  of  variance,  both  parametric  and  non- parametric, 
one-way  and  two-way.    On  the  basis  of  this  investigation,  one-way 
analysis  of  variance  is  more  useful  than  the  two-way  analysis, 
largely  because  all  the  data  can  be  used.    There  is  no  necessity  for 
an  equal  number  of  readings  in  each  column,  and  no  need  to  sub- 
stitute possibly  misleading  values  for  missing  data.    This  also 
eliminates  the  need  for  deciding  which  one  of  several  readings  in 
a  particular  month  to  use.    However,  two-way  analysis  is  useful 
in  determining  the  nature  of  the  row  effects.    There  were  cases  where 
significant  differences  occurred  between  rows,  which  indicated  a 
significant  difference  from  year  to  year,  though  this  was  rare. 
These  effects  are  shown  in  Table  3B  .  When  two-way  analysis  is 
used,  the  highest  values  per  month  should  be  used  except  in  the 
case  DO  where  lowest  values  should  be  used.    If  a  value  is  missing 
the  mean  value  should  be  substituted. 

One  may  be  tempted  to  use  non- parametric  tests  as  the  main 
method  of  analysis  because  of  several  seeming  advantages  including 
ease  of  computation  and  the  assumption  of  no  particular  frequency 
distribution.    However,  if  the  data  is  from  a  well-behaved  distri- 
bution, we  throw  away  information  by  using  the  non- parametric  tests 
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Instead  of  parametric  tests.    This  may  mean  we  can  prove  a  signifi- 
cant difference  between  two  or  more  groups  by  a  parametric  technique 
but  not  by  a  non- para  metric  technique.    However,  it  is  still  best  to 
perform  some  non- parametric  tests  along  with  the  parametric  tech- 
niques. 

In  case  studies  like  this  more  Information,  possibly  obtained 
using  recording- type  meters,  is  required  to  determine  the  population 
distribution  of  the  variables.    Once  this  is  known,  the  validity  of  the 
various  forms  of  testing  can  be  implicitly  determined.    This  kind  of 
population  distribution  information  is  required  for  all  the  parameters 
to  provide  the  best  possible  understanding  of  the  analytical  techniques, 
and  may  be  worthwhile  obtaining  sometime  in  the  future. 

Based  on  the  short-term  data  used  in  this  study,  the  following 

methods  are  best  applied  to  the  analysis  of  the  data: 

Parametric:  Analysis  of  Variance      —      one-way,  two  or  more 

unequal  samples,  using 
all  availabe  data. 
two-way,  two  or  more 
equal  samples,  using 
high  values  (except  in 
case  of  DO). 

Non- parametric:  Mann- Whitney  U-test  —     one-way,  two  unequal 

samples,  using  all 
available  data. 
Kruskal-Wallls  Analysis  of  Variance 

one-way,  three  or  more 
unequal  samples. 
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In  this  study,  all  the  tests  mentioned  above  gave  generally  consistent 

results,  even  though  the  nature  of  the  tests  varied  to  some  extent. 

Though  the  grouping  of  months  varied  from  one  parameter  to 

another,  definite  patterns  were  discernible.    This  case  study  revealed 

that  it  was  possible  to  produce  a  grouping  of  months  which  apply  in 

general  to  all  eight  parameters: 

February,  March,  April 
May,  June 
July,  August 
September,  October 
November,  December,  January. 

The  above  grouping  could  be  applied  to  other  streams  and  stations  to 
confirm  the  results  of  this  case  study.    Further  study  is  under  consider- 
ation to  support  the  validity  of  this  grouping. 

This  study  showed  that  it  is  possible  to  reduce  the  number  of 
monthly  readings  required  for  the  particular  station  studied  thereby 
achieving  some  financial  gain.    At  the  same  time  no  information  which 
is  required  for  decision- making  is  lost.    Of  course,  consideration 
should  also  be  given  to  social,  economic  and  political  factors  involved 
in  making  the  final  decision  on  the  nature  of  the  monitoring  program. 
Based  on  the  statistical  analysis  of  this  short-term  data,  it  is  possible 
to  reduce  the  number  of  samples  required.    Of  course  one  must  remember 
that  any  investigation  based  on  short-term  data  is  not  as  valid  as  that 
based  on  long-term  data. 
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The  analysts  presented  here  Is  based  on  temporal  readings  but 
it  can  be  applied  to  spatial  values  as  well;  that  is,  finding  significant 
differences  between  readings  at  different  places,   say,  on  a  lake.    As 
more  information  on  the  population  distribution  becomes  available, 
different  testing  methods  could  be  applied  which  may  change  the 
grouping  pattern.     However,  at  present,  it  is  possible  to  apply  the 
above  techniques  to  other  streams  and  stations  with  some  mathematical 
confidence. 


26 

APPENDIX,  *  NOTATION 

a  =       level  of  significance,  probability  of  a  Type  I 

error; 

BOD  =       5-day,  2Q°C,  Biochemical  Oxygen  Demand; 

o 

DO  =       Dissolved  Oxygen 

F  =       the  F  test:    parametric  analysis  of  variance; 

P  =       F  =    MSA     =     (m-1)  SSA     (rows); 

*  A  a   '(k-l),(n-l)(k-l)  "     MSE  SSE 


F 


MSA 


F      ,  m    MSB     =     (k-l)SSB       (columns); 


B  a     (n-l),(k-l)(n-l)         MSE  SSE 

k  =       the  number  of  rows; 


mean  sum  of  squares  (rows)      *    SSA  ; 

k-1 


SSB 
MSB  =       mean  sum  of  squares  (columns)  -  -!^ —  ; 

n-  1 

MSE  =       mean  sum  of  squares  (residual  or  error)  = 


SSE 


(n-l)(k-l) 

fj,  =  the  population  mean; 

n  =  the  number  of  samples  (columns); 

N  =  the  total  number  of  random  variables; 

NQ3+N02  =  nitrate  plus  nitrite; 

NSD  ~  no  significant  difference; 

ppm  =  parts  per  million; 

SD  =  significant  difference; 

SSA  ■  sum  of  squares  (rows); 
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SSB  ■  sum  of  squares  (columns); 

SSE  =  sum  of  squares  (residual  or  error); 

SST  =  sum  of  squares  (total); 

T  =  Student's  t  test:    a  parametric  test; 

x,  x..  ■  the  random  variables; 

x  =  the  sample  mean. 
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G  R  0  U  P  I  N  G 


MATER 

QUALITY 
WWAMnTERS 

RATIO: 

RES 

U  L  T  S 

GROUPS 

RANGE 

HIGH 

PARAMETRIC 

NON 

-  PARAME1 

RIC 

LOk 

1  -  MAY 

KRUSKAL- 
WALLIS 

mm  - 

WHITNEY 

MALD  - 
HOLFOUITZ 

DISSOLVtD 
UXVGEN 

FEB.,  MAR.,  APR. 

5.0,  14.0 

2.80 

N  S  D 

N  S  D 

* 

* 

MAY,  JUN. 

8.0,    13.0 

1.62 

N  S  D 

* 

N  S  D 

NSD 

JUL.,  AUG.,  SE^T. 

7.0,   12.4 

1.77 

N  S  D 

N  S  D 

* 

* 

OCT.,  NOV., 

4.0,   13.0 

3.25 

N  S  0 

* 

N  S  D 

S  D 

DEC,  JAN. 

6.0,  16.0 

2.CT 

N  S  D 

• 

USD 

NSD 

COL  1  FOR* 

FEB.,  MAR,,  APR. 

812,  21000 

99.10 

H  S  D 

N  S  D 

* 

* 

MAY,  JUN. 

ICO,     7000 

70.00 

N  S  D 

* 

NSD 

NSD 

JUL.,  AUG. 

170,   14000 

82.40 

N  S  D 

• 

N  S  D 

NSD 

SEPT.,  OCT. 

252,  33000 

143.00 

N   SD 

* 

NSD 

NSD 

NOV.,     DEC,  JAN. 

280,     8000 

28.60 

N  S  0 

NSD 

* 

* 

TUTAL  SOLIDS 

FtB,,  WW.,  APR., 

270,      968 

3.56 

N  S  D 

N  SD 

* 

» 

MAY,  JUN. 

234,      400 

1.71 

N  SD 

* 

NSD 

S  D 

JUL.,  WUfi. 

264,      726 

2.75 

NS  D 

* 

NSD 

NSD 

SEPT.,  OCT. 

222,      306 

1.79 

N  S  D 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

290,      494 

1.71 

USD 

H  S  D 

* 

* 

BOD  (5  DAY) 

FEB.  MAR.,  APR. 

0.9,  22.0 

24.50 

N5D 

N  S  D 

• 

* 

MAY,  JUN.,  JUL., AUG. 

1.0,  32.0 

32.00 

N  S  0 

N  S  D 

* 

* 

SEPT.,  OCT. 

0.2,     3.2 

16.00 

N  S  D 

* 

NSD 

NSD 

NOV.,  DcC,  JAN. 

0.4,     3.3 

9.00 

N  S  D 

N  S  D 

« 

* 

SUSPtNDED 
SOLIDS 

FEB.,  MAR.,  APR. 

10,      753 

75.30 

MSB 

N  S  D 

* 

* 

MAY,  JUN. 

5,         44 

8.a 

N  S  D 

• 

NSD 

NSD 

JUL.,  AUG. 

13,         44 

3*39 

N  S  D 

* 

NSD 

NSD 

SEHT.,  OCT. 

5,         50 

10.00 

S  D 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

5,       130 

26.00 

N  S  D 

N  S  D 

• 

* 

. 

I 
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TABLE  JL.  (CONTINUED) 
STATISTICAL       RESULTS        WITH         GROUPING 

FOR       ALL  VALUES 
DUFF  IN        CREEK 


WATtR 

QUALITY 

PARAMETERS 

GROUPS 

RANGE 

RATIO: 

HIGH 
LOW 

R  t  S  U  L  T  S 

PARAMETRIC 

NON 

-  PARAMETRIC 

1   -  WAY 

KRUSKAL- 

MANN  - 

WALD  . 

WALLIS 

WHITNEY 

WiLFUWlTZ 

NITRATE  & 

- 

NITRITE 

FEB.,  MAR.,  APh. 

0.12,  5.00 

41.60 

NSD 

NSD 

* 

* 

KAY,  JUN..JUL,  AUG. 

0.00,  0.47 

- 

N  SD 

NSD 

* 

• 

SEPT.,  OCT. 

0.00,12.50 

- 

NSD 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

0.13,   1.28 

9.85 

NSD 

NSD 

* 

* 

KJELDAHL 

NITROGEN 

FEB.,  MAR.,  APR., 

0.32,  2.00 

6.25 

NSD 

NSD 

* 

m 

MAY,  JUN.,  JUl.,  AUG. 

0.00,    1.70 

- 

NSD 

NSD 

* 

» 

SEPT.,  OCT. 

0.26,  0.78 

3.00 

NSD 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

0.30,  0.78 

2.80 

NSD 

N  S  0 

* 

* 

TOTAL 

PHuSi-HATE 

FEB.,  MAR.,  APR. 

0.08,  2.10 

26.20 

NSD 

NSD 

* 

* 

MAY,  JUN. 

0.07,  0.90 

12.90 

NSD 

* 

NSD 

NSD 

JUL., AUG. .SEPT., OCT. 

0.04,19.0 

475.00 

NSD 

NSD 

* 

• 

NOV.,  DEC,  JAN. 

0.02,  0.58 

26.00 

NSD 

NSD 

* 

* 

NSD  -        NO      SIGNIFICANT     DIFFERENCE 

S  D  -        SIGNIFICANT      DIFFERENCE 

*  -        TEST  NOT  APPLICABLE. 
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TABLE 


STATISTICAL        RESULTS        WITH  GROUPING 

(USING     MAXIMUM     VALUES) t 
DUFFIN        CREEK 


WATER 

QUALITY 

PARAMETERS 

R     E     ! 

>    U    L     T 

S 

. 

P  A  B 

a  M  E  T  R   1  C 

RON* 

P  A  R  A  M  E 

T  R    1   C 

GROUPS 

l-WAY 

2  - 

WAY 

FRIEDMAN 
2-  WAY 

KRUSKAL  - 

NALLtS 

1    -  WAY 

MANN  - 
WHITNEY 

WALD   - 
WOLFOr.lTZ 

DfFF.TJr 
COL.  B's 

&IFF.    IN 
ROWS     A'< 

DISSOLVED 
OXYGEN 

FfcB.,  MAR.,  APR. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

* 

MAY,  JUN. 

N  S  D 

S  D 

S  D 

* 

* 

NSD 

NSD 

JUL.,  AUG.,  SEPT. 

NSD 

NSD 

NSD 

NSD 

NSD 

<* 

* 

OCT.,,  NOV. 

N  S  D 

NSD 

N  S  I 

• 

* 

NSD 

SD 

DEC,  JAN. 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

COL  1  FORMS 

FEB.,  MAR.,  APR. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

* 

MAY,  JUN. 

NSD 

NSD 

NSD 

• 

* 

NSD 

NSD 

JUL.,  AUG. 

NSD 

NSD 

NSD 

• 

* 

NSD 

NSD 

SEPT.,  OCT. 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

NSD 

NSD 

N  SO 

NSD 

NSD 

* 

* 

TOT  L 
SOLIDS 

FEB.,  MAR.,  APR. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

* 

MAY,  JUN. 

NSD 

NSD 

NSD 

* 

* 

NSD 

S  D 

JUL.,  AUG. 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

SEPT.,  OCT. 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

NOV.,  DEC,  JAN. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

* 

5-DAY 

FEB.,  MAR.,  APR. 

NSD 

NSD 

NSD 

NSD 

NSD 

• 

* 

HAY,  JUN.,  JUL..AUG. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

• 

SEP ! . ,  OC1. 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

i 

NuV.,  DEC,  JAN. 

NSD 

NSD 

S  D 

N  S  D 

NSD 

* 

* 

SUSPENDED 
SOLIDS 

FEB.,  MAR.,  APR. 

NSD 

NSD 

NSD 

NSD 

NSD 

• 

• 

MAY,  JUN. 

NSD 

NSD 

NSD 

• 

* 

NSD 

NSD 

JUL,     AUG 

NSD 

NSD 

NSD 

* 

* 

NSD 

NSD 

SEPT.,  XT. 

SD 

NSD 

NSD 

* 

* 

NSD 

NSD 

NIV.,  DIC,  JAM. 

NSD 

NSD 

NSD 

NSD 

NSD 

* 

* 

i  EXCttT      IN     CAst     OF     D  0     WHERE      LOW     VALUES     WEhE     UStO. 
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TAB  It.      t  A  (continued) 

STATISTICAL         RESULTS         WITH        GROUPING 
(USING     MAXIMUM     VALUES^ 

y  u  h    f     I   II         (    k  I    t    K 


MATER 
QUALITY 

PARAMETERS 


NITRATt  & 
NITRITE 


KJELDAHL 
NITROGEN 


TOTAL 

PHOSPHATE 


GROUPS 


FEB.,  HhR.,  APR. 
MAY,  JUN.,  JUL,  AUG. 
SErT.,  OCT. 
NOV.,  DEC,  JAN. 

FL8.,  MAR.,  APR. 
MAY,  JUN.,  JUL.,  AUG. 
SEPT.,  OCT. 
NOV.,  DEC,  JAN. 

FEB.,  MAR.,  APR. 
MAY,  JUN. 

JUL.,  AUG.,  SEPT.,  OCT. 
NOV.,  DEC,  JAN. 


RESULTS 


PARAMETR    IC 


I -WAY 


DIFF.   IN 

coUiul 


N  S  D 

N  S  D 
USD 
NSD 

N  S  D 
N  S  D 
N  S  D 

NSD 

N  S  D 
N  S  D 
N  S  D 
N  SO 


2  -  WAY 


IFF. IN 


NSD 
NSD 
NSD 
NSD 

NSD 
NSD 
NSD 
NSD 

NSD 
NSD 
NSD 
NSD 


Ll .. 


*  S  D 
S  D 

NSD 
S  D 

NSD 
NSD 
NSD 
NSD 

NSD 
NSD 

NSD 
NSD 


N  0  N 


A  R  A  M  E    I   R   I  C 


FK I EOMANN 


NSD 

NSD 

» 

NSD 

NSD 
NSD 

* 

NSD 

NSD 

* 

NSD 

NSD 


KRUSHAL- 
WALUS 


NSD 

NSD 

» 

NSD 

NSD 

NSD 

* 

NSD 
NSD 

NSD 

NSD 


MANN  • 
WHITNEY 


NSD 

* 


NSD 

* 


NSD 

* 


W0LFWIT7 


NSD 

* 


NSD 

* 


NSD 

* 


NSD 

S  D 

» 


if 


NO      SIGNIFICANT     DIFFERENCE 

SIGNIFICANT  DIFFERENCE 

TEST  NOT    APPLICABLE 

EXCEPT    IN   CASE    OF  D  0     WHERE   LOW  VALUES  fcERE  USED. 


TABLE 


A  N  t 

L  Y  5 

1  S 

0  F 

VARIANCE         (I 

-HAY) 

(USING    ALL     VALUES) 

a  -  o.05 

• 

DUFF   IN        C 

REEK 

WTER  g^LlTY 

NO.     OF 

SUM     OF 

DEGREES 

"^7 

MEAN 

FA 

F 

READINGS 

OF 

PARAMETER 

GROUPS 

PER     MONTH 

SQUARES 

FREEDOM 

SQUfcRE 

CAL. 

TAB. 

REMARKS 

DISSOLVED  OXYGEN 

FEB.,  M*R.,  APR. 

8 

7 

7 

4.22 

2 

2.11 

0.324 

3-52 

N  S  D 

MAY,  J UN. 

6 

6 

4.94 

II 

4.94 

2.215 

4.96 

N  S  D 

JUL.,  AU6.,   S£t-T. 

9 

7 

6 

3.34 

2 

1.67 

I.I  10 

3.52 

N  SD 

OCT.,  NOV. 

7 

6 

1.44 

1 

1.44 

0.242 

4.84 

X  S  D 

DEC,  JAN. 

6 

7 

19.40 

1 

19.40 

2.665 

4.84 

N  S  D 

COL  1  FORMS 

FEB.,  MAR.,  APR. 

8 

7 

7 

8 
1.14  x   10 

2 

5.71  x  10 

3.23} 

3.52 

N  S  D 

MAY,  JUN. 

5 

5 

6.56  x  106 

I 

6.55  x  10 

1.644 

5.32 

N  SD 

JUL.    ,  AUG. 

9 

7 

2.C2  x  I07 

1 

2.02  x  I07 

0.859 

4.60 

N  S  D 

SEPT.,  OCT. 

6 

7 

8 
1.34  x   10 

1 

8 
1.34  X   10 

1.509 

4.84 

N  S  D 

NOV.,  DEC,  JAN'. 

5 

6 

7 

6 
7.54  x  10 

2 

3.77  x  10 

0.589 

3.68 

N  SO 

TOTAL  SOLIDS 

FEB.,  MAR.,  APR. 

8 

7 

7 

3.54  x  I04 

2 

4 
1.77  x  10 

0.898 

3.5S 

N  S  D 

MAY,  JUN. 

6 

5 

15.7 

1 

15.7 

0.008 

5.12 

N  SD 

JUL.,  AUG. 

9 

7 

6.56  x   I03 

1 

6.58  x  10 

0.489 

4.60 

N  SD 

SkPT.,  OCT. 

6 

6 

616 

1 

616 

0.313 

4.96 

N  S  D 

NOV.,  DEC,  JAN. 

6 

6 

7 

3.92  x   I03 

2 

1.49  X   10 

0.523 

3.63 

N  SD 

BOD  (5-OAY) 

FEB.,  MAk.,  Apr. 

8 

7 

7 

34.3 

2 

34.3 

0.833 

3.52 

N  S  D 

MAY,  JUN.,  JUL,  AUG. 

6 

5 

9 

7 

92.8 

3 

30.9 

0.904 

3.03 

MSB 

SEPT.,  OCT. 

6 

7 

UC4 

1 

1.04 

1.228 

4.84 

N  SD 

NOV.,  DEC,  JAN. 

6 

6 

7 

0.246 

2 

0.123 

0.126 

3.63 

N  S  D 

m 

TABLE  2     'vC°*TIVUEQ) 

ANALYSIS        OF         VARIANCE  (    I     -    W  A  Y  ) 

(USIMG    ALL     VALUES) 

DUFF  IN         CREEK 


taATER  QUALITY 
PARAMtFE'. 


SUSPEND  d) 
SOLIDS 


NITRATE    4 

NITRITE 


KJELDAHL 
NITROGEN 


TOTAL 
PHOSPHATE 


GR'JUP 


No.     OF 

READING 
Per    month 


FEB.,  KmR.,  APR. 
MAY,  JUN. 
JUL.,  AUG. 
SEPT.,  OCT. 
NOV.,  DEC,  JAN. 

FEB.,  MAR.,  APR. 
M«Y,  JUN.,  JUL.,  AUG. 
SEPT.,  OCT. 
NOV.,  DEC,  JAN. 

FEB.,  MAB.,  APR. 
MAY,  JUN.,  JUL.,  AUG. 
SEPT.,  OCT. 
NOV.,  DEC,  JA?J. 

FEB.,  MAR.,  >R. 
MAY,  JUN. 

JUL.,  AUG.,   SEPT.,  OCT 
NOV.,  DEC,  JAN. 


6 
6 

9 
8 

6 

B 

6 
B 

6 

8 

5 
5 
6 

e 

5 

s 

5 


SUM     OF 
SQUARES 


9.99  X   10 
I09.O 
89.9 
637.0 
1.39  x  10 

1.996 
0.090 
12.2 
0.194 

0.371 
0.159 
0.0734 
0.0634 

0.814 
0,0172 
39.400 
0,072 


DEGREES 

OF 
FREEDOM 


2 
I 
I 
I 
t 

2 
3 
I 

2 

2 
3 

I 
2 

2 
I 

3 
2 


MEAN 
SQUARE 


4.34  X  10 

109.0 

89.9 

697.0 

693.0 

0.998 
0.010 
12.200 


CAL. 


0,0971       0.765 


0.135 
0.0531 
0.0734 
0.0267 

0.407 
0.0172 
13.100 
0.038 


1.960 
0.982 
0.768 
5.816 
0.740 

1.046 
0.505 
1.034 


1.296 
0.534 
3.803 
1.259 

2.236 
0.163 
1.084 
1.136 


TAB. 


3.52 

4.90 
4.84 
3.74 

3.52 
3.03 
4.84 
3.63 

3.52 
3.05 
4.96 
3.83 

3.59 

5.59 
2*99 
3.74 


R>..WilG 


N  S  D 
N  S  D 
N  S  D 
S  D 
N  S  D 

N  S  D 
NSO 
N  S  D 
N  S  D 

N:S  D 
N  S  D 
NSO 
N  S  D 

N  S  D 

N  S  D 

N  S  D 

N  S  D 
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TABLE 


2  A 


STATISTICAL 


PARAMETERS         (USING     MAXIMUM     VALUES) 
DUFF    IN         CRELK. 


WATER 

qw-l  i  ry 

PAhAMETtRS 

GROUPS 

RANGE 

KAIIO: 

HIGH 
LOW 

SUM 

0  F 

SQUARES 

DEGREE 

0  F 
FREEDOM 

VARIANCl      OF     ANALYSIS 
ONE  -  WAY 

"MEAN     SQUARE 

F 
CALC  1 

F 

TAB 

MSB 

M  S  E 

DISSOLVED 

■ 

OXYGEN 

FEB.,  MAX.,  APR. 

5.0,    13.8 

2.76 

4.14 

2 

2.07 

3.90 

0.53 

3.89 

MY,   JUN., 

8.0,    13.0 

1.63 

8.80 

1 

8.80 

1.69 

5.21 

5.32 

JUL.,  AUG.,  SEPT. 

7.0,   10.6 

1.52 

1.50 

2 

0.75 

1.62 

0.46 

4.26 

OCT.,     NOV. 

4.C,    I3.C 

3.28 

3.86 

1 

3.86 

0.33 

4.63 

5.99 

DEC,     JAN. 

6.0,    13.2 

2*30 

15.13 

1 

15.13 

2.56 

5.30 

5.32 

COL  1 FOKHS 

FEB.,  MAR.,  APR., 

230,  21000 

91.3 

8.43x1 0' 

2 

4<22xt 

2.34x  id 

1.80 

3.39 

MAY,     JUN. 

100,     7000 

70.0 

7J29XICI 

1 

a77xl^ 

avTuio 

1.93 

5,32 

JUL.,  AUG., 

170,   14000 

82.1 

2.79x  1  0? 

1 

ZTB%\( 

2.97x10 

0.94 

5.99 

SEPT.,  OCT. 

260,  36000 

138.0  . 

2.49x1  O8 

1 

l.26xl( 

a4)xin 

1.96 

5.99 

NOV.,  DiC,  JAN. 

420,     8000 

19.0 

7 

1.62x10 

C, 

B.10xK 

6.52x10 

1.24 

4.26 

TOTAL  SOLIDS 

FEB.,  MAR.,     APR, 

306,       962 

3.14 

9.39  xlO4 

2 

4.95xlC 

2.07x10 

2.38 

3.89 

MAY,     JUN. 

234,      400 

1.71 

0.50 

1 

0.5 

2.08x10 

0.0002 

5.32 

JUL.,  AUG. 

274,      726 

2.65 

5.83XI03 

1 

5832 

2.40x16* 

0.24 

5.99 

SEPT.,     OCT. 

292,       398 

1.36 

682.87 

1 

633 

1625 

0.42 

7.71 

NOV.,  DEC,     JAN. 

342,      494 

1.45 

1 ,92x  Id3 

2 

958 

3055 

0.3t 

4.26 

BOO  (5  DAY) 

FEB.,  MAR.,     APR. 

O.S,     22.0 

24.40 

29.28 

2 

14.64 

28.70 

0.51 

3.39 

W*Y,  JUN.,  JUL., AUG. 

1.0,     32.0 

32.00 

132.60 

2 

44.20 

53.50 

0.76 

3.49 

SEPT.,  OCT. 

0,9,       3.2 

3.56 

0.03 

1 

0.03 

0.56 

0.05 

5.99 

NOV.,  DEC,     JAN. 

2.1,       3.6 

1.71 

0.21 

1 

0.11 

0.21 

0.50 

4.26 

SUSPENDED 

1.55x10 

i*    3.39 

SOLIDS 

FtB.,  MAR.,     APR. 

19,       753 

39.6 

2 

77.49 

29.41 

2.94x1  ( 

MAY,     JUN. 

10,         44 

4.4C 

240 

1 

240.05 

74.20 

3.23 

5.32 

JUL.,     AUG. 

24,         44 

1.83 

91.13 

1 

BUI 

39.96 

2.28 

5.99 

SEPT.,  OCT. 

13,        50 

3.85 

394 

1 

394.33 

112.28 

3.51 

\}m\fO' 

NOV.,     DEC,     JAK. 

16,       130 

8.12 

2.05x1(5 

2 

LCExii: 

IJBxri 

0.80 

5.41 

NITRATE  & 

NITRITE 

FEB.,     MAR.,  APR. 

0.51,  5.00 

12.50 

1.54 

2 

0.77 

1.54 

0.59 

3.89 

MAY,  JUN.,  JUL.,  AUG, 

0.00,  0.47 

- 

0.005 

3 

0.002 

0.029 

0.06  . 

3.49 

SEPT.,     OCT. 

0.00,12.50 

- 

18.  17 

1 

18.75 

18.17 

0.97 

5.99 

NOV.,  DEC,  JAN, 

0.46,   1.28 

2.79 

0.06 

2 

0.03 

0.10 

0.27 

4.26 
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TABLE       2  A  (con'inu£d) 
STATISTICAL         PARAMETERS         (USING  MAXIMUM  VALUES) 

DUFF   IN         C  R  E  F.   K 


MATER 

QUALITY 
PAkAMETtRS 


KJELDAHL 
NITROGEN 


TuTAL 
PHOSPHATE 


GROUPS 


FtB.,  MAR.,  APR. 
MAY,  Jon.,  JUL,,  AUG. 
SEPT.,  OCT. 
NOV.,  DEC,  JAN. 


FEB.,  MAR.,  APR. 
MAY,  JUN. 

JUL.,AUG.,SEHT.,OCT. 
NOV.,  DEC,  JAN. 


RANGE 


0.40,  2.00 
0.00,  1.70 
0.38,  0.78 
0.30,  0.78 

0.16,  2.10 
0.07.0.SO 
0.11,19.0 
0.12,  0,68 


RAT1U: 
HIGH  } 
LOW     I 


SUM 


T 


DEGREE 
0  F   LME*K  SQUARE 


r 


5.00 

2.17 

2.60 

13.10 
12.9 
173.0 
4.83 


SOWRESi  FREEDOM    M  S  Bj  K  S  E  j  CALC   ' 


TAB. 


0.18 

0.45 
0.07 
0.08 

0.98 

0.02 

66.85 

O.lb 


2 
3 
I 
2 


3 

2 


0.09 
0.15 
0.07 
0.12 

0.49 

0.02 

21.95 
0.08 


0.16 
0.13 
0.01 
0.03 

0.20 

0.16 

21.83 

0.01 


0.56 
1.21 

5.01 
4.51 

2.46 
0.10 
1.00 
5.34 


3.89 
3.49 
5.99 
4.26 

3.39 
7.71 
3.49 
4,26 


TABLE 


3A 


VARIATION         BE 

T  W  E  E  N 

COLUMNS 

(B'S) 

TWO 

WAY        ANALYSI 

S 

a  • 

■    0.05 

DUFF  IN 

C  R  i 

.  E  K 

WATER     QUALITY 
PARAMETER 

GROUPS 

No.  UF 
READINGS 
PER  MONTH 

SUM    OF 

SjUA:<ES 

(sse  ) 

DEGRtES             MEAN 

OF               SQUARE 
FREEDOM             (MSB) 

SUM     OF 
SQUARE 
(SSE) 

F     CALC. 

F     TAB. 

REMARKS. 
D IFF. BETWEEN 
COLUMNS  fl'V 

TWO-*AY 

TWO-WAY 

DISSOLVED 

OXYGEN 

FEB.,  MAR.,  APR. 

5 

4.35 

2 

2.17 

42J3 

0.475 

4.46 

N  S  D 

MAY,  JUN. 

5 

3.84 

1 

3.84 

1.26 

11,174 

7.71 

S  D 

JUL.,  AUG.,  SEPT. 

4 

2.75 

2 

1.37 

10.99 

2.835 

5.14 

N  SD 

OCT.,  NOV. 

4 

1.13  X  10 

1 

1.13  X  I0"2 

27.60 

0.003 

10.10 

N  SO 

BEC,  JAN. 

5 

4.50 

1 

4.50 

156.8 

0.458 

7.71 

N  S  D 

COL  1  FORMS 

FEB.  MAR.,  AeR. 

5 

1.05  x   10 

2 

5.25  x  I07 

8 
1.96  x  10 

2.197 

4.46 

N  S  D 

MAY,  JUN. 

5 

6 
8.45  X   10 

1 

8.45  x   10 

1.19  x  10 

2.841 

7.71 

N  S  D 

JUL.,  AUG. 

4 

2.79  x   I07 

1 

7 
2.79  X   10 

3,37  x   10 

2.482 

10.10 

N  S  D 

SEkT.,  OCT. 

4 

2.22  x   10 

1 

2.22  x   I08 

4.08  X   I08 

1.634 

10.10 

N  S  D 

NOV.,  DE-.   JAN. 

4 

7 
1.37  X   10 

2 

6.87  x  10 

2.05  x  10 

1.993 

5,14 

HSU 

TOTAL  SOLIDS 

FEB.,  MAR.  APR. 

5 

9.30  x   I04 

2 

4.65  x  id* 

1.53  x  10 

2.435 

4.46 

N  S  D 

MAY,  JUN. 

5 

t,70 

1 

1.70 

6.91  X   103 

0.001 

7.71 

N  S  D 

JUL.,  AUu. 

4 

5.83  X   10 

1 

5.83  x   I03 

8.55  x   I04 

0.204 

10.10 

N  S  D 

SEPT.,  OCT. 

3 

683 

1 

683 

5.40  X  JO 

0.253 

18.50 

N  SO 

NOV.,  DEC,  JAN. 

4 

3 
2.90  x   10 

2 

1.45  x   10 

4 

I. 41  x   10 

0.817 

5.14 

NSD 

BOD   (5  -  day) 

FEB.,  MAR.,  APR. 

5 

44.10 

2 

22 

172.60 

1.022 

4.46 

N  S  D 

MAY,  JUN.,  JUL.,  AUG. 

4 

131 

3 

43.7 

493.83 

0.797 

3.86 

N  S  D 

SEPT.,  OCT. 

4 

0.245 

' 

0.245 

3.28 

0.224 

10.10 

N  S  D 

MOV.,  DEO,  JAN. 

4 

-2 
7.17  X   10 

2 

3.58  x   lO"2 

0.63 

0.339 

5.14 

N  S  D 

TABLE      3  A     (continued) 

VARIATION    6ETI.EEN    COLUMNS    (B»s) 

a  -  0.05 
DUFF  IN    CREEK 


TWO 


WAY 


ANALYS  IS 


WATER     QUALITY 

No.  OF 

SUM     OF 

DEGREES 

MEAN 

SUM    OF 

F    CALC. 

F     TAB. 

REMARKS 

READINGS 

OF 

SQUARE 

BIFF,  between; 

PARAMETER 

GROUPS 

PER  MONTH 

SQUARES 

FREEDOM 

SQUARE 

(SSE)     , 

TWCMAY 

TtaO-»AY 

SUSPENDED 

1.62  x    I06 

SOLIDS 

FEB.,  MAR.,  APR. 

5 

2 

8.09  X   I04 

2.42  X   I05 

2.878 

4.46 

N  SD 

may,  JUN. 

5 

292 

1 

292 

30? 

3.210 

7.71 

MSB 

JUL.,  AUG. 

4 

91.1 

1 

91.1 

137 

1.990 

10.10 

N  S  D 

SEPT.,  OCT. 

4 

406 

1 

406 

349 

3.487 

10.10 

N  SD 

NOV.,  DEC,  JAN. 

3 

2.05  x   I03 

8 

(.02  x   I03 

2.90  x   I03 

1.413 

6.94 

N  S  D 

NITRATE    4 

NITRITE 

FEB.,  MAR.,  APR. 

5 

3.14 

2 

1.57 

11.3 

1.199 

4.46 

N  S  D 

MAY,  JUN.,  JUL.,  AUG. 

4 

9.13  X  I0"3 

3 

-3 

3.04  x   10 

0.1 1 

0.258 

3.86 

N  S  D 

SEPT.,  OCT. 

4 

21.3 

1 

21.3 

87 

1.120 

10.10 

NSO 

NOV.,  DEC,  JAN. 

4 

0.195 

2 

9.74  x   10 

0.31 

3.094 

5.14 

M  SD 

IUELDAHL 

NITROGEN 

FEB.,  MAR.,  APR. 

5 

0.343 

2 

0.172 

1,60 

0.860 

4.46 

N  S  D 

MAY,  JUN.,  JUL.,  AUG. 

4 

0.496 

3 

0.165 

1.23 

1.211 

3.86 

N  S  D 

SEPT.,  OCT. 

4 

6.48  X   I0"2 

1 

6.48  x   I0*2 

0.06 

4.042 

10.10 

N  S  D 

NOV.,  DEC,  JA*. 

4 

3.07  x   10 

2 

1.54  x  ID*"2 

0.09 

1.081 

5.14 

N  S  D 

TOTAL 

PHOSPHATE 

FEB.,  MAR.,  -PR. 

5 

1.13 

2 

0.565 

1.48 

3.080 

4.46 

USD 

MAY,  JUN. 

3 

1.60  x   I0"8 

1 

I.6C  x  I0"2 

0.48 

0.007 

18.50 

N  S  D 

JUL,,  AUG.,   SEfT.,  OCT. 

4 

67.5 
7.43  X  I0-2 

3 

22.5 

-2 
3.72  x   10 

197.90 

1.023 

3.86 

N  S  D 

NOV.,  DEC,  JAN. 

4 

2 

O.tl 

2.051 

5.14 

USD 

TAB  L  E 


3  B 


VARIATION 


BETWEEN    ROWS   (A's  )  -    TWO-WAY 

a-  0.05 
DUFFIN    CREEK 


ANALYSIS 


W.JER  QUALITY 
PARAMETER 


GROUPS 


No.  OF 

READINGS 
PER    Gr.CUP 


I 


SUM     OF 


SQUARES 


DEGkElS 

OF 
FREEDOM 


MEAN 

SQUARE 

(M  S  A) 


MEAN 
SQUARE 
(M  S  £) 


F     CALC. 


TWO-v<AY 


F     TAB. 


TWO-WAY 


REMARKS 
01 FF. BETWEEN 
ROWS  A's 


DISSOLVED 
OXYGEN 


COL  I  FORMS 


TOTAL  SOLIDS 


BOD   <5-DAY) 


FEB.,  MaR.,  APK. 

MAY,  JUN. 

JUL,,  AUG.,  SErT. 

OCT.,  NOV. 

DEC,  JAN. 

FEB.,  MAR.,  APR. 

MAY,  JUN. 

JUL.,  AUG. 

SEPT.,  OCT. 

NOV.,  DEC,  JAN. 

FEB.,  MAR.,  AhR. 

MAY,  JUN,, 

JUL,,  AUG. 

SEPT.,  OCT. 

NOV.,  DtC,  JAN. 

FEB,,  WAR.,  A^R. 

MAY,  JUN.,  JUL.,  AUG, 

SEPT.,  OCT. 

NOV.,  DEC,  JAI'i. 


15 
10 
IS 

8 

e 

lb 
10 

a 

8 
12 
IS 

IC 
8 
6 

IS 

l£ 

16 
6 
12 


11.03 

12.6 

3.7 

37.80 

28.10 

9. 10  X  I07 

I .84  x   I07 


8 


I  .45  x  10 
8.78  *  I08 
4.07  x  I07 
1.07  x  I05 


9.74  X   10 
5.85  x  10 


4 


I. 10  x  10 
I .36  x  I04 

173 

208 
0.625 
1.760 


4 

4 
3 
3 

4 
4 
4 
3 
3 
3 
4 
4 
3 
2 
3 
4 
3 
3 
3 


2.5 

3.17 

1.2 
12.6 

9.35 
2.27  x  I0? 
4.59  x  I01 


e 


4.82  x   10 


6 


1.26  x  10 
1.36  X  I07 
2.68  X  I04 
3.44  x  10 

1.95  x  10 

548 
ft.63  x  I03 

43.2 
69.5 

0.208 
0.537 


5.27 
0.31 
1.33 
4.35 
9.83 
2.39  X  10 
2.97  x  I06 

1.12  x  10 

8 


6 


! .36  x  10 

3.42  x  10* 

1.91  x  I04 

1.73  x  10 

2.86  x  IC4 

2.70  x  I03 


2.36  x   10 

21.57 

54.87 

1.09 

0.  II 


1.747 

5.27 

11.901 

6.39 

3.092 

4.76 

2*899 

9.28 

0.951 

6.39 

0.952 

3,94 

1.543 

6.39 

4.290 

9.28 

0.928 

9.28 

3.961 

4.78 

1.402 

3.34 

1.41 1 

6.39 

0.632 

9.28 

0.203 

19.00 

1.923 

4.76 

2.005 

3.84 

1.263 

3,86 

0.190 

9.28 

5.551 

4.76 

N  S  D 

S  D 
USD 
USD 
N  SD 
N  SD 
N  SD 
N  S  D 
N  S  D 
N  S  D 
N  S  D 
USD 
N  S  D 
N  S  0 
N  S  D 
N  S  D 
N  S  D 
N  S  D 

S  D 


to 


T  A  6  L  E  3.8  (continued) 


VARIATION 


BETWEEN     ROWS    (A's)    -    TWO-WAY    ANALYSIS 
DUFFIK    CREEK 


|    WATER     QUmLITY 

NO.     OF     1 

SUM     OF 

DEGRELS 

MEAN 

MEAN 

F    CALC. 

F     TAB. 

REMARKS 

PARAMETER 

GROUPS 

RE/JJINGS 
PER  GROUP 

SQUAKES 

■ 

OF 

FREEDOM 

SQUARE 
(MSA) 

SQUARE 

(m  S  e) 

DIFF.BtTWEEN 
ROWS  A'S 

TWO-WAY 

TWO-WAY 

SUSPENDED 

SOLIDS 

FEB.,  MAR.,  APR. 

15 

1.21   x   I05 

4 

3.02  X  10 

3.04  X  I04 

1.006 

3.84 

MSB 

MAY,  JUN. 

10 

247 

4 

61.6 

90.85 

0.679 

6.39 

N  S  D 

JUL.,  AUG. 

8 

102 

3 

34.1 

45.79 

0.745 

9.28 

N  S  D 

SEPT.,  OCT. 

8 

334 

3 

11 1.0 

116.46 

0.957 

9.28 

N  S  D 

NOV.,  DEC,  JAN. 

9 

4.78  X  10 

2 

2.39  x  10 

724.24 

3.302 

6.94 

N  S  D 

NITRATE  & 

NITRITE 

FEB.,  MAR.,  APR. 

15 

3.64D 

4 

0.910 

1.40 

0.696 

3.34 

H  S  D 

HAY,, JUN.,  JUL.,  AUG. 

16 

0.242 

3 

8.C6  x  10*^ 

0.01 

6.835 

3.86 

S  D 

SEPT.,  OCT. 

8 

56.200 

3 

18.5 

190.00 

0.974 

9.28 

N  SO 

NOV.,  DEC,  JAN, 

12 

0.504 

3 

0.168 

0.05 

5.343 

4.76 

SD 

KJELDAHL 

NITROGEN 

FEB.,  MAR.,  APR. 

15 

0.321 

4 

8.C2  x  KT 

0.20 

0.401 

3.64 

N  S  D 

HAY,  JUN.,  JUL.,  AUG. 

16 

0.263 

3 

8.78  x   10** 

0.14 

0.644 

3.30 

N  SD 

SEPT.,  OCT. 

8 

4.11   x   I0"2 

3 

1.37  x  I0"2 

0.02 

0.854 

9.28 

NSD 

NOV.,  DEC,  JAN. 

12 

0.159 

3 

5.27  x  \6* 

0.01 

3.712 

4.76 

N  S  D 

TOTAL 

PHOSPHATE 

FEB.,  MAR.  APR. 

15 

0.953 

4 

0.238 

0.18 

1.290 

3.34 

NSD 

MAY,  JUN. 

6 

0.162 

2 

8.12  X   lO-2 

10.24 

0.341 

19.00 

NSD 

JUL.,  AUG.,  SEPT.,  OCT. 

16 

64.000 

3 

21.3 

22.00 

0.97C 

3.86 

NSD 

NOV.,  DEC,  JAN. 

12 

0.117 

3 

0.389 

0.02 

2.145 

4.76 

NSD 

o 


FINAL  GROUPING 
FIGURE  I  -  DUFFIN  CREEK 


Significant  Difference  (2— way  parametric,  high  values) 
Significant  Difference  (1-way  parametric,  all  values) 


